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Abstract 

The operation of a source of entangled electron spins, based on a superconductor and two quantum 
dots in paralleii, is described in detail with the help of quantum master equations. These are derived 
including the main parasitic processes in a fully consistent and non-perturbative way, starting 
from a microscopic Hamiltonian. The average current is calculated, including the contribution of 
entangled and non-entangled pairs. The constraints on the operation of the device are illustrated 
by a calculation of the various charge state probabilities. 
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I. INTRODUCTION 



Entanglement is a basic resource in quantum computation and quantum communication^. 
Recently, various experiments for quantum information processing schemes have been suc- 
cessfully implemented with photons as Bell inequality violation^ or teleportation^. 

Any system with a two-level quantum degree of freedom is a possible candidate to carry 
a quantum bit. One of such is the electron and its spin. In principle, individual electrons 
can be manipulated in a quantum circuit and have the advantage of promising high-level 
integration in electronic devices^. Notice that the electron flow can be in principle much 
larger than the photon flow in equivalent optical devices where attenuation is necessary to 
produce individual photons. Moreover, photons essentially do not interact except during 
their generation process, whereas Coulomb correlations between electrons in a quantum 
circuit open the possibility for new operations between quantum bits^. 

Non-locality in quantum mechanics can be probed by letting two quantum degrees of 
freedom interact, and subsequently separating these two systems. Here, electronic entangle- 
ment can be created using a superconductop^*^, where two electrons forming a Cooper pair 
are in a singlet state. The superconductor is coupled to two arms, each of them collecting 
one electron from each Cooper pair. The emission of one electron in each lead from the 
same Cooper pair corresponds to the so-called Crossed Andreev proces a^°i^^i^^ , which can 
be understood as a non-local Andreev reflection: the emission of one of the electrons can 
be seen as the absorption of a hole with opposite spin and opposite momentum. The two 
electrons forming the singlet are then spatially separated. It is then necessary to avoid the 
"ordinary" Andreev reflection where the two electrons go into the same lead. This selection 
can be enforced, either with the help of spin fllters, leading to energy entanglement^. Or 
alternatively, one can use energy filters, leading to spin entanglemeniii^. Quantum dots 
with Coulomb blockade, inserted in each branch, can efficiently select the crossed Andreev 
process. As another possibility, the superconductor can be replaced by a normal quantum 
doiji^. In this paper, the studied device consisting of a superconductor connected to two 
quantum dots in parallel will be called the Entangler (see FigHj). Branching currents in 
the right and left leads were calculated for this Entangler in Ref. Q using a T-matrix ap- 
proach. Entanglement can be probed by sending the electrons from a splitted pair into a 
beamsplittepi^ and by measuring noise correlations^^. 
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In the present paper a microscopic derivation of quantum master equations^^ for the 
Entangler is presented. It provides a simple, intuitive approach to probe entanglement 
and to monitor the effect of parasitic processes. Compared to a T-matrix derivationi, this 
approach has the advantage of describing the whole charge dynamics in a non-perturbative 
way (this statement will be qualified below). This allows to derive not only the average 
current but also the higher moments of the current distribution. Another point is that 
quantum master equations can be applied to any arbitrary quantum system containing 
superconducting elements, or to another kind of Entangler. 

Over the past years a great interest has been devoted to the description of the trans- 
port properties through devices containing coupled nanostructures, where quantum inter- 
ference has a strong influence. A rather accessible method, generalizing the classical master 
equationsii, has been developed in Refill where Bloch-type quantum rate equations have 
been derived using the Schrodinger equation. When the system is an isolated quantum dot 
in the Coulomb blockade, only the diagonal elements of the density matrix (the occupation 
probabilities) enter the rate equations. On the other hand, when the transfer of electrons 
through a quantum device goes through a superposition of states in the different parts of this 
device, non-diagonal matrix elements will appear in the equations of motion. The master 
equations then take into account coherent processes and are a generalization of the Bloch 
equationsi^. 

The microscopic derivation of these equations provides a good understanding of the cor- 
respondence between quantum and classical descriptions of transport in mesoscopic systems. 
The crucial point is the decoupling between the time scales which specify, first, the dynamics 
inside the reservoirs and, secondly, the inverse rates for coupling the quantum states and the 
leads. This decoupling procedure is justified as long as the time scales characterizing transfer 
within the quantum system and injection (emission) from (to) the reservoirs are both large 
compared to the time scale for fluctuations within the reservoirs. This is equivalent to a 
markoffian hypothesis^^. 

Until now, the generation of quantum master equations has been limited to the case of 
sequential tunneling within quantum dots coupled to normal reservoirs, using a microscopic 
Hubbard-type Hamiltonia n^^i^" . In particular, electron transfer which is mediated by high 
lying virtual states is not accounted for. Consider now the case of a superconductor coupled 
to quantum dots: single-electron tunneling does not conserve energy and is forbidden as 
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the electron transfer is accompanied by the emission of a Bogolubov quasiparticle. However 
two-electron events such as Andreev processes (transfer of a pair of electrons out of the 
superconductor) and superconducting cotunneling (S-cotunneling) processes^! (transfer of 
an electron from one dot to another via the superconductor) connect low energy states, and 
thus enter the lowest-order contribution to the tunneling current from the superconductor. 
One simplification would be to assume that the two-electron tunneling processes occur si- 
multaneously, and are described by a pair Hamiltonian: rate equations have been written 
recently in this manner for the transport processes in a teleportation cell which employs 
an array of normal and superconducting quantum dots^. However, in presence of transport 
channels mixing different processes, it is safer to derive quantum master equations starting 
directly from the microscopic Hamiltonian. This is achieved in the present work, taking 
into account the main parasitic processes. The sequence of relevant steps will clearly require 
virtual states which contribute to Andreev and cotunneling events. After having established 
the equations including coherent quantum mechanical effects and Coulomb blockade, we will 
determine their range of validity and show the relevance of the lifetime of virtual states. The 
derivation of quantum master equation is non-perturbative with regard to transitions within 
the Entangler, while the coupling to the leads is treated within the Fermi golden rule as in 
the orthodox theory of Coulomb blockade^^. 

The paper is organized as follows. In SeclUl we present the system and its energy scales, 
together with the Crossed Andreev process - the main process - and the important parasitic 
processes that can occur during its evolution. This allows to write the many-excitation 
wavefunction which is the starting point of each microscopic derivation. This derivation is 
first described in Sec lHIl for the Crossed Andreev process, without any parasitic process. 
Parasitic processes are presented next, and compared in Sec IIVI before giving the complete 
description of the system by quantum master equations in SecEl and appendix [B1 Sec lVIl 
provides the physical discussion of the operation of the device as a function of its parameters. 
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II. THE ENTANGLER DEVICE AND ITS PARAMETERS 



A. The model 

Let us first provide a qualitative description of the Entangler. The setup involves a 
superconductor (S) coupled by tunneling barriers to two quantum dots {Di and D2) which 
are themselves coupled to normal leads, L and R (see FigC]). Only one level is retained in 
each dot, assuming the level separation in each dot to be large enoughi. The energy levels of 
the dots can be tuned by external gate voltages. The microscopic Hamiltonian of the entire 
system is the following: 

TC = TLq + TCtunnel (1) 

where 

TYo = ^ Ek'jl^^ka+Eid\^dia+E2dl^d2a+Uini„ni^+U2n2an2^+y~^^ Eia]^ai„+^ Er-al^ara 

k,a la ra 

(2) 

where ^ka, di^, ai^, a^a are destruction operators for Bogolubov quasiparticles, dot electrons 
and reservoir electrons, = d\^di^ is the occupation number in the dots, which enters the 
Hubbard repulsion term with coupling constants Ui and U2. A possible inter-dot repulsion is 
omitted here for sake of simplicity, but it could easily incorporated in the energies of various 
charge states of the two dots system. 

The tunnel Hamiltonian which connects these elements by a one-electron transition reads: 



'Htunnel = ^ ^kldl^Cka + ^ ^~k2dl_^C^k~a + ^ ^laja^la + ^ fir-aL'^^T + h.C. (3) 

k,a k,a l,a r,a 

with a single electron tunneling amplitude Qi {Q2) between S and Di {S and D2), and 
fi; (Qr) between Di and L (between D2 and R). a = {|,— |} is the spin variable. Note 
that Htunnel is writtcu in the Fourier space. Point contacts are assumed between 5* and 
dots 1 and 2 (in fi and thus the tunneling term is Qidj^Crc^, which can be written in 
the Fourier space: J^k^i^^''" Aa'^ka = J2k^kidl„Cka- The effective momentum dependence 
of the tunneling amplitude 17^, introduces a geometrical factor, which can strongly influence 
the transition amplitude for processes involving the two quantum dots. During the injection 
process, Cooper pairs are initially separated into one electron in a dot and one quasiparticle 



5 



in 5*. We introduce the Bogulubov transformation: 
with 




FIG. 1: The Entangler setup: a superconductor injects electrons in quantum dots Di and Z?2, whose energies 
in state |1) (i.e. one excess electron) are respectively Ei and E2. Electrons in the dots can subsequently 
tunnel into the normal reservoirs L, R. 

Here S stands for the annihilation of a Cooper paip2^, and (ps is the superconductor's 
phase. The two electrons from a Cooper pair become an entangled pair of electrons (only the 
singlet state is involved) when going into different leads. Current flow is imposed by a voltage 
bias A/i between the superconductor and the leads. The basic mechanism for entanglement is 
based on a Crossed Andreev process between the superconductor and the two quantum dots, 
forced by the Coulomb blockade in the dots. First, two entangled electrons are created in 
Di and D2 via a virtual state which contains a quasiparticle in S whose energy is larger than 
As, the superconducting gap. This process is coherent, and couples the superconducting 
chemical potential ns and the final energy of the pair in the dots Ei + E2. This Andreev 
process probability is optimized at Ei + E2 = fis, and behaves like a narrow two-particle 
Breit-Wigner resonance. Then the two electrons tunnel independently to each lead. This 
whole sequence of events forms the Crossed Andreev channel. 
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B. Working conditions 

Next, the relevant parameters describing the device are discussed, following Refill. First, 
the charge states of the quantum system have to be well separated to avoid transitions due 
to thermal excitations. But the thermal energy must be large enough in comparison to 
the transition probability to allow the markoffian hypothesis. Therefore Tl^r -C ksQ -C 
Ei — fJ'L,R- In order to conserve spin and thus the singlet state during the electron transfer, 
spin-flip must be excluded. Thus each dot cannot carry a magnetic moment which could 
interfere with an electron coming from S, i.e. it must carry an even number of electrons^. 
Moreover, when an electron is deposited on a dot, another electron of this dot with opposite 
spin could in principle escape to the normal leads thus spoiling the entanglement. This 
spin-flip process is suppressed when the dot level spacing 6e is larger than the imposed bias 
A/i and the temperature A;^©. Entanglement loss can also occur because of electron-hole 
excitations out of the Fermi sea of the leads during the tunneling sequence. Such many- 
particle contributions are suppressed if the resonance width ^ = 2'itpljj{Ei2)\^l,r{Ei^2)\'^ 
is smaller than £'12 — fJ'L,R- This justifies the microscopic Hamiltonian of EqlH 

Next, given this Hamiltonian, one needs to justify the derivation of the quantum master 
equation. Single-electron tunneling from the superconductor to the leads via the dots is 
avoided because it implies the creation of a quasiparticle in S. This process costs at least 
As which is assumed to be much larger than Afi and ksQ. 

C. Parasitic processes 

The main purpose of this device is to force the two electrons from a pair to propagate 
in the two different leads. In a clean three-dimensional superconductor, this process is 
decreased by a geometrical factor 7^ = 6 '"^o ^^^j:'^'^^^ (^0 is the superconductor coherence 
length and r = |ri — is the distance of the two contacts between dots and S). The 
crossed Andreev amplitude is then 'JaT, with T = {7i/2)N{0)QiQ2- Beside the decay on 
.^0, the algebraic factor can be improved by reducing the dimensionality^^ or using a dirty 
superconductop24. Incidentally, the finite width of the contacts may introduce diffraction 
corrections to the geometrical factor. Note that when taking into account the finite thickness 
of the contacts, the geometrical factor can be modified^^. 
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There are three main parasitic processes which could decrease the Entangler efficiency. 
Two of them create different channels of emission of two electrons coming from a Cooper 
pair, for which the two electrons can tunnel to the same lead^. Although they involve 
higher energy intermediate states, those do not suffer from the geometrical factor of the 
Crossed Andreev channel. In addition, an elastic cotunneling - this process will be called 
S-cotunneling in what follows - connects every channel to other processes by transferring 
an electron between the two dots via S. 

The two electrons of a Cooper pair can tunnel through the same dot by an Andreev process 
(cf. FigOI). Because of double occupancy, the pair would get an energy U due to Coulomb 
repulsion. This is a coherent process between two energy levels with a large energy difference 
U. Because this energy cost is much larger than the Andreev process probability amplitude 
Tj ~ N{0)Qf involving a single lead, this process is strongly suppressed^. Alternatively, a 
pair could propagate to the same lead if the first electron injected on a given dot leaves it 
before the second electron is deposited on either dots. It goes to the corresponding lead 
while its twin electron "has been staying in S"' as part of a quasiparticle (cf. FigEj). The 
latter can then choose toward which dot it will tunnel. It will prefer the same dot in order to 
avoid paying the geometrical factor. This latter process costs Ag and thus can be suppressed 
with As ^ 'JaT. Let us notice that this process requires three transitions, including one 
transition to a reservoir, thus it is not coherent. 

By a S-cotunneling process via S, an electron can tunnel from Di {D2) to D2 (-Di)(cf. 
Fig. E}. This is a coherent process between two discrete energy levels, Ei and E2 for a 
single electron in the two dots or Ui and Ei + E2 for a doubly occupied dot (f/2 and Ei + E2 
for the opposite configuration). Cotunneling is characterized by an amplitude 7c 7", with its 
own geometrical factor 7c . If the energy difference between the two coupled levels is much 
larger than the process amplitude T^, this process will be weak. 

To summarize, the working regime of the device is the following: 

A5, U, 1^1 -E2\>5e> A/i, kBQ > Tl^r, Ta, Tc (8) 

This working regime contains the justifications for the approximations made in the derivation 
of the master equation: the markoffian approximation and the relevant processes involving 
at most two successive virtual states with only one quasiparticle in S. 
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In what follows we also assume that A/i > ksQ, in order to ensure the irreversibility of 
the pair production. 



III. MASTER EQUATIONS FOR THE CROSSED ANDREEV CHANNEL 



The transport channels which are described above can be characterized by the charge 
configuration of the isolated quantum system for each step of the Entangler operation. The 
quantum system is composed of the dots and the superconductor, but its dynamics can be 
directly probed by integrating out excitations in the reservoirs and superconductor. Using 
the Schrodinger equation and generalizing the procedure of reflisl. it is shown here how to 
derive quantum master equations which describe the evolution of the reduced density matrix 
of the system. As a starting point, we consider the dynamics in the situation where only 
the Crossed Andreev process and one-electron relaxation processes are effective - the ideal 
regime. The wavefunction is thus chosen to include only the charge states involved in this 
particular channel. A reduced Hilbert space containing the lowest energy states and the 
required virtual intermediate states is chosen (containing a single quasiparticle in S). 

state V b|i^ state c h^rSg 



state a 



dot 1 



dot 2 



o 
o 



Left 



vRight 



0.2 



a 



O 



O2 



Oi 



state b 



Or 



Ol 



Oi 



state v' b2ka 



state d b2isc5 



Or 



^ state a bij-s^. 



FIG. 2: Sequence of states for the Crossed Andreev channel of the Entangler. For instance, bik,a- denotes 
the amplitude to have an electron in dot i while a quasiparticle is created in the superconductor. First an 
electron is deposited in either dot, next the second electron tunnels and forms a singlet state in the pair of 
dots, next either electron is absorbed in the reservoir, and finally the two dots are empty. 
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The many-excitation wavefunction for this problem is written as: 



+ J2lr,a blrS,a{t)al^al^S 



|0) 



(9) 



+ 



+ 



+ 



where h,,,{t) are the time-dependent amphtudes for finding the system in the corresponding 
states with the initial conditions feo(O) = 1 and all other 6(0) are zero. The indices indicate 
the electron occupation in the dots and reservoirs, as depicted in Fig. El The use of 
Schrodinger equation and the form of call for some comments. In fact, as said above, 

the temperature is not zero thus one should in principle rely on a density matrix description 
from the beginning. Yet, under the condition F^ /j < fc^G < A/i, one can simply use the 
Schrodinger equation in a reduced subspace of states^^. Those states for instance do not 
include electron-hole excitations in the same reservoir: these are supposed to relax on a very 
short time, due to inelastic processes occurring in L and R. On the contrary, all possible 
charge and spin states on the dots, together with all excitations including holes in L and 
electrons in i?, are considered. Summing on these reservoir states eventually lead to the 
equations for the two reduced density matrij^i^. 

After substituting Eq.ljHI) into the Schrodinger equation i\ip{t)) = Ti.\ip(t)), an infinite 
set of coupled linear differential equations is obtained for b{t) by projecting i{ipi\ip(t)) = 
{■ipi\T-C\ip{t)) for each state \ipi). \ipi) characterizes the quantum state of the total system 
including the environment. Applying the Laplace transform 



^0 

and taking into account the initial conditions, an infinite set of algebraic equations is ob- 




(10) 
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tained for the amplitudes b{E) (see FigiS}: 

{E + ir])bo - ^ avkQlJ)ika + ^ (rvk^*_kJb2ka = i (Ha) 

ka ka 

{E + irj- El- Ek)hika = crvlVtkiho + Uk^%J)i-2Sa (Hb) 
{E + ir] - E2- Ek)b2ka = -(TvlQ-k2bo - Uk^^libusa (He) 

{E + irj -El- E2)bi2Sa = ^ U*k^-k2blka - ^ ulVtklhka + ^ ^MrSa + ^ ^^/^2iSa 

fe k r I 

(lid) 

{E + irj - El - Er)birSa = ^rbusa + ^ ^^z&ZrS^ (He) 

{E + i7]-E2- Ei)b2lSa = ^ibuSa + ^rkrSa (llf) 

r 

[E + irj -El- Er)birSa - ^ 0-'Vk^*klblrSa,lk'a' + ^ Vk^*^k2blrSa,2k' a' = ^ihra + ^rhla 

k'a' k'a' 

(llg) 

Each term corresponds to the transition between two successive states. Each transition 
leads to the creation or annihilation of a quasiparticle either in S* or in a reservoir. There 
is an fundamental difference between the two types of transitions. The first one involve 
an excited state whose lifetime is so small (r^p ~ l/A^ ^ ^/T) that coherence is kept 
until the quasiparticle is destroyed. On the other hand, in the reservoirs, quasiparticles 
instantaneously decay {rreiax ~ ^/Ep <^ l/F) so coherence is lost (Markoff process). To 
simplify the system of equations, the expression for b is substituted in terms of the type 
from equations containing sums. Every sum over the continuum states {k, I, r) is 
replaced by integrals (see AppendixEI)- Crossed terms (like '^ibiQiflr/ {E — Ei)) vanish^^, 
and the following set of equations is obtained: 

[E + irj- 2c{Ti + T2))h = 2^ATe'^%bi2s,a - bi2s,-a) (12a) 

{E + i7j-Ei-E2- c\Ti + T2) + £f)bi2S.a = 2a^ATe-'Ho (12b) 

[E + lTj-Ei- Er + tY)blrS,a = ^rbl2S,a (12c) 
{E + lTj- E2- Ei+ l^)b2lS,a = ^lbl2S,a (l2d) 
(E + irj - El - Er - 2c{Ti + T2))blrS,a = 2'^ATe^'^^ {blrSa,l2S'a> - blrSa,l2S'-a'^ + ^lblrS,a + ^rb2lS,a 

(12e) 
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with Ti = jN{0)Qj and c, c' are numerical constants (see Appendix EI), involved in self- 
energy corrections. Here the coefficients for virtual states (states \v) and in FiglSI) have 
disappeared from the equations. This is the consequence of the succession of quasiparticle 
creation and annihilation transitions forced by the assumption that two quasiparticles cannot 
coexist in S. 

The singlet/triplet basis is now chosen. For instance, in the global wave function, 
h2sAA-. is replaced by b-f^\d\A^-d\A2a)/V2 + b'^^^^^^^ 
From Eq. (|12b|) one can say that coefficients &i2S,cr and &i25,-cr for a given spin are op- 
posite. This is the same for birs,a and birs-a, ^2«s,(t and &2/s,-a- The tunnel Hamil- 
tonian conserves spin, therefore there is no coupling towards triplet spin states. Thus 

7sinnlet /Ki /Ki j Ttriplet r\ 

bijs = y^bijs,a = -V2bijs,^ and 6^^/ = 0. 

The density matrix elements of the set-up are now introduced. The Fock space of the 
quantum dots consists of four possible charge states: \a) - levels Ei and E2 are empty, \b) - 
levels El and E2 are occupied, |c) - level Ei is occupied, \d) - level E2 is occupied. Reservoirs 
states are identified by n, the number of pairs of electrons out from S to the reservoirs. To 
obtain the reduced density matrix, elements are summed over n: 



E 

n=0 



a 



in) 



(13) 



In every state, electrons are paired in a singlet state. The matrix elements are defined 



as: 



an 



(^bb 



O'dd 



+ 



Lr 



^singlet 



+ 



El singlet 
\'^l'r'S',12S 

l',r' 
2 



l<l',r<r' 
2 

+ ■■• 



E 7 singlet . \ ^ j singlet 
rlrS + 2-^ ri'r'S',lrS 



l',r'<r 



I 1 2 I 

E ^singlet _|_ \ ^single 



et 
21S 



l'<r,r' 



+ 



+ 



boK^f* + ^ 



Tsingleti singlet* , 
/ , "irS "lrS,12S' + ' ' ' 



Lr 



C^ba 



ab 
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The matrix density elements are directly related to the coefficients b{E) by a Laplace 
transform: 

l...,r... 

where a/ (3 specify the charging states associated with the amplitudes (fe's). The equations 
for 77, = can be obtained straightforwardly. For instance, to get a^J , Eq. ()12a| is multiplied 
by b^lE') and the conjugate equation written for E' is subtracted. 



bb 



= 2x/2z7^T (e^^-ai? - e^^^ai")) (15a) 
(r, + r^) ai? - 2V2t^AT (e-^<^«ai? - e^^^a^^) (15b) 



= -r.a(? + Tnaf^ (15c) 
-S = -r.aS+r.aS) (15d) 
= 2V2Z7AT (e-^^a« - e^^«a«) + T^a^ + T^a^") (15e) 



Note that the diagonal matrix elements (the "populations") are coupled with the off-diagonal 
density-matrix elements ("coherences"), which is symptomatic of a coherent, reversible tran- 
sition. 

To obtain the equations for the coherence one subtracts Eq. ()12ap for E multiplied by 
k2s,singiet{E') and Eq.^ for E' multiplied by bl{E): 

cr^J = -li^L + Tn) ^i? +z{E, + E, + K{T, + T,)) a^^ + 2v^^7^Te^^^ (a^ - ^iJ)) (16) 

where K = c' — 2c. 

These equations describe the sequential evolution of the system and involve consequently 
only processes between real states. Coherent processes (not involving reservoirs) couple 
non-diagonal elements to diagonal elements while relaxation processes couple only diagonal 
elements. From the set of equations [TH one can see that these processes do not interfere 
because of the loss of phase coherence introduced by the markoffian approximation, i.e. the 
sum over reservoir states. A density matrix element for one particular state is then only 
coupled to the elements for adjacent states in the sequence. Thus the processes can be added 
easily, which will be crucial when considering the full operation including all channels. 

Here, because only one current channel is implied in the ideal operation, we can easily 
verify that the equations are the same for each n. Therefore the sum over n is obvious 
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and one obtains the master equations for the evolution of the density matrix describing the 
system: 

(Taa = 2v^z7aT (e-^^V^b - e^'^'aba) + r^a^e + TuaM (17a) 
ahb = -2V2i^aT {e-"^'aab - e^'^^a^a) - {Tl + Tr) a^b (17b) 

CTcc = -'^LCTcc + rfiCTfeft (17c) 
CTdd = -^R<^dd + ^L(^bb (17d) 

c^ab = -\ {Tl + r^) aab + i{E[ + E'^) a^b + l^ilATe''^^ {a^a - a^fc) (17e) 
with E[ = Ei + KTi. 

This is the main result of this section. First, let us remark that the transition rates 
Tl,r appear only from the dots to the reservoirs, and not in the opposite direction. This 
is consistent with the assumption that ksQ is small compared to the transition energies 
between dots and reservoirs. This limitation of Gurvitz's method is not a problem here 
since the Entangler actually needs to be strongly biased to avoid decoherence effects. The 
second term of Eq. (ll7e|] expresses that two discrete energy levels are coupled by a coherent 
process involving two transitions. Note that the probability of transmission between these 
two states is maximum in the resonant case, e.g. e = E[ + E2 is zero. 

The ideal operation of the system involves only one channel for transferring a Cooper 
pair to the reservoirs: the two electrons tunnel towards different leads. Actually, using the 
normalization condition for the populations (Taa + O'bb + O'cc + O'dd = 1, equations ^) are easily 
solved for the stationary current, I = I{t 00) (d^^ = 0): 

ir/e = T^abb + T^a^ = 87ir^ + + ^''^ 

Pj^' = eVnatb + eVnadd = iT (19) 

This current is made of entangled singlet pairs. This result was obtained earlier in Refill in 
the limit 7yiT <^ T and Tl = F^j. Here the presence of the term 87^T^ in the denomination 
comes from a complete (non-perturbative) treatment of both Andreev and decay processes. 

The equality of the currents in the two branches of the device is a direct consequence of 
the Crossed Andreev process. Every electron pair crosses and goes out of the system - each 
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electron on its own side - before the next pair is injected in the dots. Those "cycles" never 
overlap in this ideal working regime. 
In the case T ^ 'JaT, one obtains: 

Jr'/e ^ (20) 



while in the case T <^ 'JaT, 



irVe = (21) 

like a single quantum dot between two leads^i. In the latter situation, the dots are almost 
always occupied, so that the resistance is dominated by the rate associated with the two 
barriers - in parallel - between dots and leads. 



IV. PARASITIC CHANNELS 



The ideal working regime is affected by parasitic processes: Andreev tunneling via a 
single dot, one-by-one tunneling or S-cotunneling. The two first ones have been separately 
computed by the T-matrix in RefQ. Their effect is to create different channels of pair current 
which decrease the efficiency of entanglement. As said before, the terms for each process can 
be added in the equations and combined before including them together in a whole system 
of quantum master equations collecting every possible processes (see partE|)- To start with, 
the different processes will be separately considered. 



A. Direct Andreev effect process against Coulomb blockade 

Let us imagine that a Cooper pair tunnels to the same quantum dot by an Andreev 
process, while generating a doubly occupied state. 

Because of Coulomb repulsion, an energy U (Eq.^) is required for having two electrons 
in a same quantum dot. If U is "large enough", such a process will have a low probability. 
With conventional dot technology, the interaction energy U ~ IK in the quantum dots can 
be controlled so that it is smaller than the superconducting gap > 2K. Therefore the 
doubly occupied energy level has no coupling to the continuum of S quasiparticles, which 
would effectively introduce a broadening. 
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FIG. 3: A current channel sending a pair of electrons to a same reservoir. Andreev process towards one 
quantum dot can happen against strong Coulomb repulsion U. 



Similarly to the case of the ideal working regime (Sec lIIIj) . the set of differential equations 
associated with this Direct Andreev channel are established for the reduced density matrix 
elements. Here, only one branch - say L - is considered for simplicity. The Fock space of the 
quantum dots consists here of three possible charge states: \a) - both dots are empty, |e) dot 
1 is doubly occupied, |c) - dot 1 is singly occupied. The wavefunction takes the following 
form: 

+ j:jn{t)dld\^S 

+ ••• + ••• + 

From the Schrodinger equation, and performing steps similar to Sec lIIH the set of equations 
for the density matrix elements is: 



|0) 



(22) 



2r^(Tee 



O-r, 



dee = -2?ri (e-^<^«(T,e - €^^^(7,^) ' ^V'^^Te 

--t[U, + K'T,] aae + 2iT^e^^ (a,, - ^ee) + r'^cr. 



(23a) 
(23b) 
(23c) 
(23d) 



with K' a numerical constant, and = 2npi{Ui + Ei)\VLi{Ui + -Ei)]^ the level broadening 
introduced by coupling of the two-electrons level with lead L. These equations are similar 
to equations (jB)). Nevertheless, the sequence passes through a high energy-level (U) via 
an Andreev process which implies an oscillation between two discrete energy levels, and 
U[ = Ui + K'Ti. On the contrary, in the ideal regime, this energy difference can be as small 
as desired. 
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B. One-by-one electron tunneling to the reservoir 

This channel is another way to send a pair into one single lead. Before the second 
electron of a broken Cooper pair can tunnel to a dot, the first one already leaves the dot to 
the corresponding lead. The second electron will tunnel through the same dot as its twin 
electron with a much higher probability (FigEI) than through the other dot, because of the 
geometrical factor. The latter process will be simply neglected. 



State a state c state a 




FIG. 4: Sequence corresponding to the tunneling of a singlet pair through one branch of the device. States 
I a) and |c) are coupled through two successive virtual states. 

There are only two processes involved in this channel. The first one, between states 
I a) and |c), requires two consecutive virtual states, both containing a quasiparticle in S. 
Because of the coupling with a continuum of states in the lead, phase coherence is lost thus 
off-diagonal matrix elements - or coherences - are not coupled to populations. Therefore 
this channel is peculiar in the sense that it is incoherent even though it involves transitions 
with S. The equations describing the evolution of the density matrix are obtained as before. 
The Schrodinger equation gives: 



{E + iri)bo = i + ^ avk^l^hka 

k,(j 

{E + irj-Ei- Ek) hka = crvinkibo + ^ Vl^kka 

I 

{E + 17]- El- Ek) blka = ^iblka + UkQlJ)lla 

{E + ir]- El - El) bii^ = ^ ulflkikka + ^ cyQiibwa 

k V 

{E + ir] - El - El,) bw^ = ^li'bua - ^ibw^ + ^ (y'vk^lJ)w^^ika' 

k,cr' 



(24a) 

(24b) 

(24c) 
(24d) 

(24e) 
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Let us eliminate bi^^. To simplify, the notation Aij = E + ir] — Ei — Ej is introduced. 



Ai; -cil^ + i—\ bii„ - 2^ — — .r^(A)x ^o + ^ — : 7^ — tTTau ^1''- (26) 



Finally, using integrals calculated in Appendix lA 5[ eq. IA14[ the following set of equations 
is obtained: 

2 

, ttAc" 

Ai; - c"nl + t^- ) hi, = a^e-^'^^ibo + cj^^nj > ] n.bu'a (28) 



E + i7]-2cnl + i{^] TL\bo = t (27) 



\ 27~' '^T'' 



Aw - 2cQl + I ( ^) 1 bi,, = Qi,b,i, - Uibu'a (29) 



ttA 



Virtual states have disappeared from the equations. The remaining term in TiVti corresponds 
to the three-step process coupling \a) to |c). Introducing the elements of the density matrix 



K -^Tiaaa + '^lcTcc (30a) 



one gets : 

^cc = -Tl(Tcc + K"-^TL(Taa (30b) 

(30c) 

where K" = A/tt"^. This process behaves as the transport through a single dot where the 
first barrier between the left lead and the dot is a three-step process via two virtual states 
and the second barrier is a classic tunnel barrier. 



C. S— cotunneling between the two quantum dots 

Another process involves intermediate virtual states of the quantum device which are 
common to the other processes: cotunneling^i between the two quantum dots via S. This 
process involves oscillations between two position states and connects all of the channels 
studied until now. 
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FIG. 5: Cotunneling between the two dots. An electron from dot 1 tunnels towards the dot 2 via a virtual 
intermediate state containing a quasiparticle. Two contributions participate to the cotunneling depending 
on when the initial electron is transferred. 

It can occur in different situations: between states containing only one or two electrons 
in the two dots. Like for the Crossed Andreev process, the transmission probability depends 
on the energy difference between the two coupled states. The equation of evolution for the 
density matrix describing oscillations between two states - |c) (electron on dot 1) and \d) 
(electron on dot 2) are established: 



(Tdd = ilcT {(Tdc - (Ted) 

&cd = iAEacd + ijcT (dec - (Tdd) 



(31a) 
(31b) 
(31c) 
(31d) 



where AE = E2 — Ei, 'jc 



7% cosjkpr) 



kpr 



is the geometrical factor corresponding to this 



cotunneling process^^. Note that when the distance r which separates the two tunneling 
locations is zero, 7c diverges. This is expected because this process has no meaning for 
the same tunneling location: this local process brings back the system in the same state, it 
only participates to the renormalization of the energy level of the state by coupling with the 
continuum of quasiparticles in S. Note that the transition amplitude T is the same as for 
Andreev process. 
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V. ENTANGLER IN THE PRESENCE OF PARASITIC PROCESSES 



One of the advantage of Bloch-type equations is to be able to study all processes together 
and non-perturbatively. In the previous sections, a specific system of dynamical equations 
was obtained separately for different channels of pair current. In particular, such channels 
are repeated cycle after cycle, which allows to systematically group the contributions with 
different reservoir variables (by recurrence over the number of pairs transmitted to the leads) . 




FIG. 6: General operation including the three Andreev channels and S-cotunneHng. States with three or 
four electron states are omitted for clarity. Real states are fully squared while virtual states are dashed 
squared. To make it simpler, spin is not represented. The il's correspond to transitions between two 
quantum states \S) \dots) (g) \l,r), while Tc indicates the resonant cotunneHng process. Certain mixing 
process, like Direct Andreev effect between states |c) and state with one electron in dot 1 and two electrons 
in dot 2, are not presented for lack of space. However such processes are included in the quantum master 
equations. 

In reality, each channel (induced by Crossed- Andreev, Direct Andreev, S-cotunneling) 
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mixes into one another, so one needs to gather all transitions in a single set of equations 
for the density matrix. Because of this mixing, it is no more possible to establish a set of 
equations cycle after cycle. 

A starting point for deriving generalized quantum master equations is thus to label the 
amplitude associated with each process by variables which count how many entangled pair 
have passed through reservoir Rot L ot both (while being split). Note that such variables do 
not appear in the quantum master equation of each channel because they have been summed 
over. It is straightforward, but tedious, to write a full Schrodinger equation for the most 
general operation, combining all states, and to derive the density matrix equations. The 
basic assumption is that not more than one quasiparticle is excited in the superconductor 
during the processes. 

As was said in part IIIH all the processes can be gathered without appealing to the full 
derivation of the Schrodinger equation, by adding terms corresponding to each process. We 
set the equations for a given state of reservoirs were rii {ur) singlet pairs of electrons have 
tunneled to the reservoir L (R) and uq pairs whose electrons have tunneled to different 
leads. States are thus defined by the charge of quantum dots, one electron already in 
the reservoir while its twin electron is still in the quantum system, and ra^, and tt-q: 
\^|:) = \dotl,dot2) ® \U'rr') ® In^, n^j, rio). To get equations for only the charge states of the 
dots, they are summed over I, /', r, r' and the recurrence is made over n^, ur, uq. The 
obtained set does not depend on the number of channels and leads. Thus Gurvitz's method 
for generating quantum master equationsi^ can be generalized to the multi-terminal case 
with many current channels. The full system is given in Appendix El One can notice that 
the parasitic processes may generate triplet pairs in the leads L, R. 

VI. DISCUSSION 

The set of quantum master equations will now be used to describe more quantitatively 
the transport properties. To assess the constraints on parameters, each channel will first be 
studied, before using the complete set to obtain a numerical evaluation of the operation in 
a realistic regime. 

By solving quantum master equations one can find the average current for each uncoupled 
channel. This will be done for the symmetric case {Tl = Tr, Ui = U2 and Ti = T2) and 
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assuming that Tl^h = ^'l r = ^ when the coupling between quantum dots and lead depends 
weakly on the energy: The Direct-Andreev current is computed in the stationary regime 
with Eqs.(l23ll: 



With U' > 2T, r we have, as in Ref| 



1 fiT^ 

jAndreev ^ ^-^ __ (33) 



while with 7^T^ ^ Eq. fTKll can be written: 



jCAndreev ^ p ° M-^ /q/i\ 

~ 87iT2 + r74 ^^^^ 
The current created by the one-by-one tunneling process is given by: 

Without taking here into account elastic cotunneling, one can see here the relationship 
between parameters that must be fulfilled to approach the ideal working of the Entangler: 
U',As ^ ma,x[T,T/jA,£/lA]- This can be understood with a dynamical study of each 
channel. Actually Andreev processes are coherent processes which create an oscillation 
between the state where the Cooper pair is in S and states where the pair of electrons is in 
the dots. Thus it will be a competition between the period and the amplitude of oscillations 
and the probability of tunneling from a dot to a reservoir. Let first consider the case where 
•jaT ^ r. Then for resonant Crossed- Andreev process 

a,,{t) = 1(1- ^) (1 - cos(27ATt)) 6"^* (36) 



2 V 2T2^ 
while for Direct-Andreev process 

^eeit) = jjj^^^, (l - cos( Vf/'2 + ATH)) e-r* (37) 

Because T is small we are here in the regime where the Crossed-Andreev channel is more 
probable than the Direct-Andreev one because many oscillations between coherent states 
can occur before a transition to a reservoir has happened. On the other hand (7^?" ^ F) 
one gets: 
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As soon as the pair has tunneled to the dots, it goes to the reservoirs. And because the Direct- 
Andreev frequency is larger than the Crossed-Andreev one {U ^ IaT), there is a small time 
interval in which Direct- Andreev is favored even though the amplitude of oscillation (and 
thus tunneling between \a) and |e)) is smaller: for a relaxation time l/F of the order of half 
a period of oscillation for Direct Andreev effect {n/U), after a time t ~ l/F, the population 
of state |e) can be much larger than population of state \b). 

The same kind of argument can be given to study the effect of S-cotunneling. As said 
before, for U and As large enough, the only parasitic effect is elastic cotunneling. Us- 
ing only this process and Crossed-Andreev process in the master equation, the efficiency 
of entanglement is calculated depending on E = \Ei — which controls S-cotunneling 
probability. We want to know the proportion of electrons from a same pair tunneling to 
different reservoirs {Pentangied) or to the same reservoir (PparasiUc) ■ Cycles of current do not 
overlap so the probability is the same for each cycle. To calculate them, we can use Bloch 
equations describing the evolution on only one cycle to get first |c) and \d) populations as 
a function of time. From state \b) the first electron tunnels for example towards the left 
reservoir. The chance for the second electron to tunnel towards the right (left) reservoir is 
TRPd{t) {TLPc{t)) assuming that pd{Q) = 1. Thus Pentangied = J^TRPd{t)dt. For F^ = F^^: 

_ T^ + E^ + 2^lT- 

^entangled — p2 _|_ ^2 _j_ 4^27^2 V^^/ 

From equation (jJOl, we can see that the condition to neglect S-cotunneling, leading to 

Pentangied ~ 1, is 7cT < max[E, F] . 

A more general study using the complete set of equations (see Appendix [B| has to be 
performed. This set of equations can be solved in the stationary regime, but the general 
solution is typically cumbersome. For the sake of readability, it is presented here taking 
into account the parasitic processes only to first order. This fixes the different energy scales, 
previously discussed above, which define the working regime of the Entangler. Here the 
asymmetry F^ 7^ F^j is kept to show the role of S-cotunneling. 

h = eTiabb + eF^cTee + 2eF'^(Tcc (41) 
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FIG. 7: Charge states populations as a function of U for A5 = 9.5i^, Ei = —E2 = 0.5/^, Vj^ ji = F'^ ^ = 
T = 0.1/v, 7A,7c 0.2. States \a), \b), |c), |e), |/) refer to FigEl State \k) refers to the triplet state 
shared between dots, and states \g) and \h) refer to three electrons states (cf. Appendix^J- Iff); 1^) and 
|fc), populations correspond to the three lowest curves. The population of states containing doubly occupied 
dots vanishes when U increases. For low values oi U {U ^ \Ei — E2\), the asymmetry is introduced by 
energy difference between states |e) (two electrons in dot 1) and |/) (two electrons in dot 2). 



eao 



'aW 



-2K"Aa,T,^ (1 - ^) 
-8A.or.-(2 + ft(l-^)+§) 



(42) 



where A = ^2 .^S^-u 2 and a^^ = A + 1 + Tr/Tl + Tl/Tr. From Eq.(jl2l) we can exhibit 
which parameters are controlling each contribution to the total current. 

To complete this discussion, the set of equations is used to describe the average popula- 
tions of each state depending on some relevant parameters. A5 is taken to be the largest en- 
ergy scale. With niobium as superconductor, one takes A5 ~ 9.25K. For a two-dimensional 
quantum dot, small enough (lOnm^), one takes \Ei\ ~ 0.5K and U ~ 9K, with T, F ~ O.IK— 
and 'J A, 7c ~ 0.2. 

On figure U\ it can be seen that the population of states containing doubly occupied 
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FIG. 8: Ratio between populations of state \b) (singlet state shared between dots 1 and 2) and of state 
|e) population(two electrons in dot 1) for As = 9.5K, Ei = -E2 = 0.5K, Tl.r = T'^ ^ = T = O.IK, 
lA,lc ~ 0.2. It indicates the ratio between Direct- Andreev channel and Crossed-Andreev channel. The 
latter is strongly favored when U increases {pb/pe{U/T — 90) = 83.3). 

dots vanishes when U increases. It is important to notice that when U ~ \Ei\, the system is 
asymmetric and the channel with U ~ \Ei — E2\ is favored because a Direct- Andreev process 
becomes resonant. At the working point {U/T = 90) Pe/Pb = 0.012. Two channels can be 
compared in calculating the ratio between two populations: on figure [HI the ratio between 
the population of state \b) and the one of state |e) indicates which of Crossed-Andreev and 
Direct-Andreev process is the most likely depending on U. Thus increasing U increases the 
efficiency of entanglement. For small f/ ~ F, the two channels become comparable because 
decays to reservoirs are much faster than Crossed-Andreev oscillations. 

A large F will allow a fast transition between dots and reservoirs. That is why increasing 
F/T will favor the most likely process which connects the superconductor to the dotsi. On 
figure ini we can see that increasing F favors the decay of a single charge state before another 
pair tunnels towards the free quantum dot. Actually, because Direct-Andreev oscillations are 
faster (frequency^ ^/E"^ + AT'^) than Crossed-Andreev oscillations (frequency^ 'jaT), even 
if their probability is smaller, the decay towards reservoirs can happen before one Crossed- 
Andreev process has been achieved. Thus increasing F/T at fixed U, Direct-Andreev process 
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FIG. 9: Ratio between populations of state \h) (one electron in dot 1, two electrons in dot 2) and of state |c) 
(one electron in dot 1) for Ag = 9.5K, ^ -£2 = 0.05A', T^^r = T^^, T = O.IK, U = IK, ja:1c ~ 0.2. 
Increasing F compared to the transition rate of Direct- Andreev and Crossed- Andreev processes allows to 
favor the decay of single charge states before another Cooper pair tunnels to the free quantum dot. For 

U = IK lOT, Ph/Pe < 1.5%. 

increase to the detriment of Crossed-Andreev one. 



4 6 

r/T 
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VII. CONCLUSION 



In this article, quantum master equations have been derived, starting from a microscopic 
Hamiltonian for the superconducting-dot Entangler. Using the Schrodinger equation tech- 
nique developed in RefllSl. the full equations describing the evolution of the reduced density 
matrix are obtained, retaining as virtual states only single particle excitations in the super- 
conductor. Considering only one level by dot, all possible processes are taken into account 
in a fully consistent and non-perturbative way: Crossed-Andreev process, responsible for en- 
tanglement, as well as Direct- Andreev and one-by-one tunneling processes, and cotunneling 
through the superconductor. The latter connects all the other processes, yet the quantum 
master equations written in Appendix |B1 take into account all processes in a coherent way. 
From them, the average current has been calculated. The conditions on the Entangler pa- 
rameters, needed for an optimal operation of the device, have been derived, and extend the 
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result of Refill. 

The power of master equations is to give access, not only to the first moment, but to all 
moments of the current distributior>2i. In a forthcoming papei^, shot noise correlations are 
computed in order to give a clear diagnosis of entanglement2*22. Another extention of Bloch 
equations is to include explicitly spin/charge relaxation or coupling to external degrees of 
freedom, in order to quantitatively study decoherence effects. 

Such a derivation of quantum master equations, including higher order process, can obvi- 
ously be generalized to a wide class of quantum systems involving discrete charge states and 
coherent/incoherent transitions. It is therefore a valuable tool for investigating nanostruc- 
tures in view of controlling quantum information based on spin/charge degrees of freedom. 
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APPENDIX A: CALCULATION OF INTEGRALS 



To obtain the evolution equation of the density matrix, it is necessary to compute some 
integrals arising from the coupling between 5* and the two dots. 



1. Crossed-Andreev effect 



The tunneling of the two electrons of a same Cooper pair to two different dots gives a 



contribution (see Eqs. Illal llg 



E — Ei — Eh 



(Al) 



The two energy levels of the dots are assumed to be close to ns- The transitions amplitudes, 
VL, depend weakly on the energy so they can be considered as constant with a phase factor 
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e . Neglecting E — Ei <^ Ek ^ As one obtains: 



(27r)3y 2EI 

^1^2 V 



with /i = and the volume. 

Because of parity the integral can be extended from — oo to oo. 

I A = r dk k (e^'^ - e-^'^) ^ ^ (A3) 

~'~ \ 2m 2m 



The four poles are: 



/ 2mA \^ ±iarctan( 



k = ±kF\ll+ [j^j e ' V«^4; = ±akFe 



ki = akpe^^ , k2 = —akpe k^ = —akpe^^, ki = akpe the contour is the positive 
half-circle for e*'^^ and the negative one for e"*^''. 



{ki - k2){ki - hi){ki - ki) {k2 - ki){k2 - k:i){k2 - k^) 



(ks - ki){ks - k2){k^ - k^) {k^ - ki){ki - k2){ki - k^) 

^ 7riiiii2_v_ , / \ 7re yjakpr cose _ -iakprcose^ / 

^ 2ir (27r)3'^Uv 2(aM' sin(2^) J ^ > 

with sin26' = 2mA/{a'^h'^kj,), sin 6' = A/2Ef. Given that a and 0086* ~ 1 (A5 < ^ir), one 
obtains: 

U = |iv(o)n,n.e-^?^!ii^ (A5) 

2 kpT 
In what follows, is noted: I a = IaT with T = |A^(0)fiifi2 and 7a = e '^^^^^^f^, the 
geometrical factor for Crossed- Andreev effect. 



2. Direct-Andreev effect 

The tunneling of the two electron of a same Cooper pair to the same dot, i, gives a 
contribution 
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From the previous calculation, one must take the limit r — in the Eq. (jA5jl . The same 
result is found when making the calculation without taking into account the phase factor 
gifcx -y^^i^jf.}^ generates the geometrical factor. The amplitude of this effect towards i-side is 
then Ti = ^N{0)nl 



3. Self-energy 

The self-energy terms are due to the coupling between a discrete state (state with zero or 
one electron in a dot) and a continuum of states (quasiparticle states in S) . They correspond 
to the renormalization of of these energy levels. They involve l^^p when the annihilation of 
an electron in S corresponds to the creation of quasiparticle, and \uk\'^ when the creation 
of an electron corresponds to the creation of quasiparticle. In Eqs. ()llal llg ), this terms 
corresponds to 



E — Ei — Eu 

k * " 



for a given z-side. 

The sum are transformed into integrals over quasiparticle energies, E^^ with a density 



of states given by N{E) = N{0)E/ E"^ — A|. For the calculation, Ei/As is noted Cj and 
E/As is noted x. 
a. Terms in 

T, f 1 i\ X 



2/ VV^^ x) ei-x'^'' 

Ti , /^/^ XX ^ 67' 



In (2(1 — eA) — Ti — , ^ ( — h arcsin(ej 



^ -:^ln2 
2 

Ir = cTi (A9) 
where c is a numerical constant. 
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b. Terms in \uk\'^ : This term never appears alone, so we just have to calculate terms 
with lufcp - Ivkl"^. 




Jr = c% (AlO) 

To avoid the logarithmic divergence we introduce a physical cutoff - the electron band width 
- to get a finite result. This does not yield a large contribution because of the logarithm: if 
the band width is 1000 times higher than the gap it only gives a factor 8fi^ where Q ^ Ei. 
Self energy terms remain small. Let us define K = c' — 2c' for the following. 



4. S-Cotunneling 

Local S-cotunneling has no meaning (tunneling of an electron between two places) so 
keeping the geometrical contribution of the integrand in this process, one gets: 

With \uk\^ - \vk\^ = ^k/Ek-. 

Ic = ^I^-X- r dkkie^^^-e-^'^) S_Zi^ 

Using once again the residue theorem one gets: 

-7rO,0„ T/ /r^ \'i 1,2 -akpr sine 

/Iiilii2 V I Z,in \ /itvpt [ A{akpr cos e +9) ( 2 ie _ -ie\ _ -i{akpr cose +9) ( 2 -i9 



^ ^ r .(a/c,.cose+e) ^ 2 ^9 _ -^9^ _ -^iak,r cos 9+9) , 2 

^ 2ir (27r)3 J {akpf sin(2^) L V ; V 

(A12) 

Witha^~l + l(f|)': 



The second term is much smaller than the first one (A5 ^ Ep). The only diff'erence with 

fill instead of ^ 



the Andreev amplitude is the instead of S-cotunneling diverges for r -^^ 0. 
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5. One-by-one electron tunneling to the reservoir 



Here the calculation is not complicated by a phase factor. The sum over k is simply 
replaced by an integral over energy. 



|2 



^lk{^lk+irL/2 

El 

N{0)QlQi / —dE 



k 

r>00 



A Ji x^\/x^ — 1 
iV(0)^ (A14) 



APPENDIX B: QUANTUM MASTER EQUATIONS FOR THE ENTANGLER 

The set of fully consistent and non-perturbative quantum master equations can be derived 
(see main part of the paper). For simplicity, the space of charge states has been restricted 
here to 0, 1, 2 or 3 electrons in the two dots. Numerical calculations have been made with 
this set of equations including states \g) (one electron in dot 2, two electrons in dot 1), \h) 
(one electron in dot 1, two electrons in dot 2) and \k) (triplet state shared between dots 1 
and 2). 



baa = +2zTi {aae - (Tea) + 2zT2 ((Ta/ - (Tfa) 

2\/2i'yAT {(Jab - (Jba) + TlCTcc + ^R(ydd -2{Tl + f i?) (Jaa 

&bb = +iV2-fcT (o-fce - (Teb) + iV^lcT {Cbf - CTfb) - 2y/2i-fAT {Cab - (Tba) 
+ |rR(Tcc + ^TL(rdd + \^'L^gg + \^'R(^hh - 2 {T L + ^ r) (^bb 



(Bl) 



(B2) 



O-cc = hcT {(Ted - (Tdc) + 2iT2 {ach - (The) + 2TL(Taa + 2T RCbb + 2r'^ CTge " (r^ + 2tR)acc (B3) 



(^dd = hcT {acd - o-dc) + 2zTi (a^g - (Jgd) + ^TRCTaa + 2TL(ybb + Sr'^o-// - {Tr + 2fL)aM (B4) 
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2zTi {aae - (Tea) - iV^JcT {<Jbe - (Teb) + ^R<7gg - 21^(7, 



(B5) 




(B6) 



99 ~ 



2iTi {adg - CTgd) - (2r'i + TR)a, 



(B7) 




(B8) 




(B9) 



dab = iiE[ + E'2) (Tab + 2v/2i7^T {(Taa " O'feb) + iV2T{(Tae + C^af) " {iTia^b + ^T'sCT/b) 

-\ {2VL + 2VR + VL + Vj^aab 

(BIO) 

CTae = i{E[ + U[^)aae + 2iTi ((Taa " (Tee) +«( y27c2"crafe-2A/27^Tcrfee-2T2a/e) - (fi+fij+r'^) (Tae 

(Bll) 

(Ta/ = i{E'2 + U22)<Jaf+2iT2 ((Taa - (T//) +i (y27cT(Ta6-2A/27AT(Tfe/ -2Tiae/) - (fL+fij+r^)(Ta/ 

(B12) 

dbe = « (t^n - E'2) abe+iV2-fcT {<Jbb - cree)+i{,2Tiaba-2y/2jATaae-y/2jcTcrfe)-{TL + Tn + T^) abe 

(B13) 

^b/ = ^ (^22 - ^1) cTbf+iV2lcT {abb - aff)+i{2T2aba-2V2'jATaaf-V2^cTaef)-iTL + Tr + T^) abj 

(B14) 

CTcd = i{E2- E[) (Ted + i7c7' ((Tec - C^dd) + 22(TiCrcg - T2Crhd) - ^ + + 21^ + 2fijj (Ted 



(B15) 
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o-ch = i {U22 + E^) Cch + 2iT2{(rcc - cxhh) - hcTcxdh - {Tl + + TR)ach (B16) 

CTdg = i {U[^ + E[) adg + 2iTi{add - cTgg) - i-fcTccg - (Xr + ^'l + ^L)o-ch (B17) 

where E' and U' are the renormalized energy levels. Tl^r = K"^Tl^rTi 2/ As corresponds 
to the one-by-one process. The superconducting phase, which do not change any result, is 
omitted here. 

APPENDIX C: DERIVATION OF QUANTUM MASTER EQUATIONS FOR THE 
ENTANGLER USING A PAIR-HAMILTONIAN 

The same set of quantum master equations could be obtained from an effective Hamil- 
tonian applied to the method developed in R.eflisl. In considering all processes, this Hamil- 
tonian can be derived from the microscopic Hamiltonian (eqlH) using a projective trans- 
formation which eliminates states with quasiparticles in the superconductor to the lowest 
order: 

+2Ti(dLdU)5 + 2T2{dldl^)S + V2^cT ^1^2. 

cr 

+ J2 ^AJi, + ^ralj2a + ^AAaS + ^ ^rd\A,S + k.C. (CI) 

l,(j r,a l,a r,a 

with ^li^r = ^/,r^i,2/^s corrcspouds to the one-by-one process. The method only requires 
the amplitude for probability of processes coupling different states of the quantum system, 
and leads to a following general system: 

&al3 = i{Ei3 - Ea)(Tap + ^ I C"a7^7/3 " ^ ^ ^ayCTyP | 

7^/3 \ 77^« / 

- 2 ( ^"^T ^ Y ^f^^^ ) 2 ^ (^7^" + ^5^/3) (^y& (C3) 

\y^a 7^/3 / ySj^afi 
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where the i7's are the coherent transition matrix elements and the F's the relaxation rates. 
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